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This paper presents a new method for reducing the number of sources of evidence to combine 
in order to reduce the complexity of the fusion processing. Such a complexity reduction is 
often required in many applications where the real-time constraint and limited computing 
resources are of prime importance. The basic idea consists in selecting, among all sources 
available, only a subset of sources of evidence to combine. The selection is based on an evi- 
dence supporting measure of similarity (ESMS) criterion which is an efficient generic tool 
for outlier sources identification and rejection. The ESMS between two sources of evidence 
can be defined using several measures of distance following different lattice structures. In 
this paper, we propose such four measures of distance for ESMS and we present in details 
the principle of Generalized Fusion Machine (GFM). Then we apply it experimentally to 
the real-time perception of the environment with a mobile robot using sonar sensors. A com- 
parative analysis of results is done and presented in the last part of this paper. 

© 2010 Elsevier Inc. All rights reserved. 



1. Introduction 

Information fusion (IF) has gained more and more interests in the scientific community since the end of the 90s with the 
development of sophisticated multi-sensor and hybrid (involving human feedbacks in the loop) systems in many fields of appli- 
cations (such as robotics, defense, security, medicine). IF appears through many scientific international conferences and work- 
shops [10]. The main theories useful for information fusion are Probability Theory [18,28] (and more recently Imprecise 
Probability Theory [41 ]), Possibility Theory [7] (based on Fuzzy Sets Theory [48]), Neutrosophic Set Theory [17] and belief func- 
tion theories (such as Dempster-Shafer Theory (DST) [32] and more recently Dezert-Smarandache Theory (DSmT) [34-36]). 

In this work, we concentrate our attention on belief function theories and specially on DSmT because of its ability to deal 
efficiently with uncertain, imprecise, conflicting, quantitative and qualitative information. Basically, in DST, a basic belief 
assignment (bba) m( ) is a mapping from the power set 2® (see Section 2.2 for details) of the frame of discernment 0 into 
[0,1] such that 
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m(0) = 0 and m(X) = 1. 

Xe2® 

In DST, 0 represents the set of exclusive and exhaustive possibilities for the solution of the problem under consideration. In 
DSmT, 0 can be a set of possible nonexclusive elements and the definition of bba is extended to the lattice structures of 
hyper-power set D® (see also Section 2.3 for a brief presentation and [6,9,36] for definitions, details and examples), and 
to super-power set S® in UFT (Unification of Fusion Theories) [33,35]. In general, m( ) is not a measure of probability, except 
in the case when its focal elements (that is, the elements which have a strictly positive mass of belief) are singletons; in such 
case, m( ) is called a Bayesian bba [32] which can be considered as a subjective probability measure. In belief function the- 
ories, the main information fusion problem consists in finding an efficient way for combining several sources of evidence 
Si,s 2 , . . .,s„ characterized by their basic belief assignments (bba’s) m 1 (-),m 2 (-),- • - , m „(■)■ These bba’s are defined on the same 
fusion space, either 2®, D®, or S® depending on the underlying model associated with the nature of the frame 0. The dif- 
ficulty in information fusion arises from the fact that the sources can be conflicting (that is, one source assigns some beliefs 
to a proposition A whereas another source assigns some beliefs to a proposition B, but A and B are known to be truly exclu- 
sive (0nB = 0)) and one needs a solution for dealing with conflicting information in the fusion. In DST, Shafer proposes 
Dempster’s rule of combination, as the fusion operator, to combine sources of evidence whereas in DSmT the recommended 
fusion operator is PCR5 (Proportional Conflict Redistribution # 5) rule of combination, see [32] and [36] for discussions and 
comparisons of these rules. PCR5 is more complex than Demspter’s rule but it offers a better ability to deal with conflicting 
information. 

However, both rules become intractable in some applications having only low computational capacities (as in some 
autonomous onboard systems by example) because their complexity increases drastically with the number n of sources 
combined or with the size of the frame 0, specially in the worst case (that is, when a strictly positive mass of belief is as- 
signed to all elements of the fusion space). To avoid this problem, one can do: (1 ) to reduce the number of sources combined 
and (2) to reduce the size of the frame 0. In this paper, we propose a solution only for reducing the number of sources com- 
bined. We are not concerned about the second aspect in our application of robot perception since in this application our 
frame 0 has only two elements representing the emptiness or occupancy states of the grid cells of the sensed map of the 
environment. To expect good performances of such a limited-resource fusion scheme, it seems natural to search and com- 
bine altogether only the sources which are coherent (which are not too conflicting) according to a given measure of similar- 
ity in terms of Dedekind’s lattice. 

Such an idea has been already investigated by several authors who have proposed some distance measures between two 
evidential sources in different fields of applications. For example, Tessem [38], in 1993, proposed the distance 
dij = maX(i ie @\BetPi(8i) - BetPj(0i)\ according to the pignistic probability transform BetP(-). In 1997, Bauer [1] introduced 
two other measures of error to take a decision based on pignistic probability distribution after approximation. In 1998, 
Zouhal and Denoeux [49] also introduced a distance based on mean square error between pignistic probability. In 1999, 
Petit-Renaud [29] defined a measure directly on the power set of 0 and proposed an error criterion between two belief 
structures based on the generalized Hausdorff distance. In 2001, Jousselme et al. [13] proposed in DST framework, a new 

distance measure d,j = 1 + 2(m i , m 2 ) between two bba’s for measuring their similarity (closeness). In 2006, 

Ristic and Smets [30,31 ] defined in the TBM (Transferable Belief Model) framework a TBM-distance between bba’s to solve 
the association of uncertain combat ID declarations. These authors also recalled Bhattacharya’s distance 

= J 1 - Y^AeFiJ^BEFjVmiWnijiB} between two bba’s. In 2006, Diaz et al. [5] proposed a new measure of similarity between 
bba’s based on Tversky’s similarity measure [40]. Note that in belief function theories, the classical measures used in Prob- 
ability Theory (for example, Kullback-Leibler’s (KL) distance [2]) couldnot be applied directly because bba’s werenot prob- 
ability measures in general. Kaburlasos et al. [16,19] also proposed several valuable similarity measure functions based on an 
inclusion measure [14,15] and a metric respectively. Dongrui Wu et al. proposed a vector similarity measure for linguistic 
approximation, which was very useful in understanding the uncertainties associated with linguistic terms [45], And then, 
he also compared deeply ranking methods, similarity measures and uncertainty measures for interval type-2 fuzzy sets [46]. 

In this paper, we develop an Evidence Support Measure of Similarity (ESMS) in a generalized fusion space according to 
different lattices structures [6,9,47] for reducing the number of sources of evidence to combine. This allows us to reduce 
drastically the computational complexity of the fusion processing for the real-time implementation purpose. As shown in 
the next sections, we propose several possible measures of distance for ESMS. The purpose of this paper is neither to select, 
nor to justify, the best measure of distance for ESMS but only to illustrate from a real-time experiment with real data the 
advantages of this new approach. 

This paper is organized as follows. In Section 2, we briefly recall the main paradigms for dealing with uncertain informa- 
tion. In Section 3, we give a general mathematical definition of ESMS between two basic belief assignments and some basic 
properties of ESMS. In Section 4, we extend and present different possible ESMS functions (distance measures) fitting with 
the different lattice structures listed in Section 2. A comparison of the performances of four possible distances is made 
through a simple example in Section 5. The simulation presented in Section 6 shows in details how ESMS filter is used within 
GFM scheme and its advantage. An application of ESMS filter in GFM for mobile robot perception with real-data (sonar sen- 
sors measurements) is presented in Section 7 to show the advantages of the approach proposed here. The conclusion is given 
in Section 8. 
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2. The main paradigms for dealing with uncertainties 

2.3. Probability Theory and Bayesian rule 



Probability Theory [18] is the original theory for dealing with uncertainty. We will not present this theory in detail since 
there exist dozens of classical books devoted to it, see for example [28]. We just recall that a random experiment is an exper- 
iment (action) whose result is uncertain before it is performed and a trial is a single performance of the random experiment. 
An outcome is the result of a trial and the sample space 0 is the set of all possible outcomes of the random experiment. An 
event is the subset of the sample space 0 to which a probability measure can be assigned. Two events A, and Aj are said 
exclusive (disjoint) if A,- n Aj = 0, Vi # j, where the empty set 0 represents the impossible event. The sure event is the sample 
space 0. Probability Theory is based on Set Theory and the measure theory on sets. The following axioms have been 
identified as the necessary and sufficient condition for probability P( ) as a measure: Axiom (1 ) (nonnegativity) 0 ^ P(A) < 1, 
Axiom (2) (unity) P(@) = 1, and Axiom (3) (finite additivity, 1 ) if A^A^. . .,A n are disjoint events, then P(A t uA 2 U . . . UA n ) = 
Events which are subsets of the sample space are put in one-to-one correspondence with propositions in belief func- 
tion theory [32] (pages 35-37). That’s why we use indifferently the same terminology for set, event or proposition in this paper. 
The probabilistic inference is (usually) carried out according to Bayesian rule: 



VB.P(B)> 0, P(A,|B) 



P(Aj n B ) 
P(B) 



P(B1A,)P(A,) 



0 ) 



where the sample space 0 has been partitioned into exhaustive and exclusive events A^,A 2 ,. . .,A n , that is, A, n Aj = 0, (i # j) 
and Ai u A 2 U ■ ■ ■ U A n = <9; P( ) is an a priori probability measure defined on 0 satisfying Kolmogorov’s axioms. In Bayesian 
formula, it is assumed that the denominator is strictly positive. A generalization of this rule has been proposed by Jeffrey 
[11,12] for working in circumstances where the parochialist assumption is not a reasonable assumption, i.e. when 
P(BjB) = 1 is a fallacy, see [12,25] for details and examples. 

Using the classical terminology adopted in belief function theories (DST or DSmT) and considering for example 0 = {A,B}, 
a discrete probability measure P() can be interpreted as a specific Bayesian belief mass m( ) such that 

m(A) + m(B) = 1. (2) 



2.2. Dempster-Shafer Theory (DST) 

In DST [32], the frame of discernment 0 of the fusion problem under consideration consists in a discrete finite set of n 
exhaustive and exclusive elementary hypotheses 0,, that is, 0 = {0i,0 2l . . .,0„}- This is called Shafer’s model of the problem. 
Such model assumes that an ultimate refinement of the problem is possible, exists and is achievable, so that elements 0,-, 

i= 1,2 n are precisely defined and identified in such a way that we are sure that they are truly exclusive and exhaustive 

(closed-world assumption). The set of all subsets of 0 is called the power set of 0 and is denoted 2®. Its cardinality is 2 |e| . 
Since 2® is closed under u and all 0 f , i = 1,2,. . ,,n are exclusive, it defines a Boolean algebra. All composite propositions built 
from elements of 0 with u operator such that: 

(1) 0,0!,... ,0 n g2®; 

(2) If A, Be 2®, then A u B e 2®; 

(3) No other elements belong to 2®, except those obtained by using rules (1) or (2). 

Shafer defines a basic belief assignment (bba), also called mass function, as a mapping m( ): 2® -> [0,1] satisfying m(0) = 0 
and the normalization condition. Typically, when 0 = {A,B} and Shafer’s model holds, in DST one works with m( ) such that 

m(A) + m(B) + m(A u B) = 1, (3) 

m(A u B) allows us to assign some beliefs to the disjunction A u B which represents the ignorance in choosing either A or B. 
From this simple example, one clearly sees the ability of DST in offering a better model for a total ignorant/vacuous source of 
information which is obtained by setting m(A u B) = 1. For Probability Theory, one would be forced to adopt the principle of 
insufficient reason (also called the principle of indifference) to take m(A) = m(B) = 1/2 as default belief mass for representing 
a total ignorant body of evidence. 

In DST framework, the combination of two belief assignments mj() and m 2 () is done with Dempster’s rule of combina- 
tion. This rule is nothing but a normalized version of the conjunctive rule in order to remove the total conflicting mass [32]. 
Dempster’s rule is defined if and only if the two sources of evidence are not fully conflicting, that is, when 

J2vv i7!i(Xi)m 2 (X 2 ) 5^ 1 and it is mathematically defined by m(0) = 0 and for X ^ 0 by 
Ai,a 2 6 z 

X, nX 2 = 0 



Another axiom related to the countable additivity can be also considered as the fourth axiom of Probability Theory. 
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J2x, x 2e 2 ®m,(X 1 )m 2 (X 2 ) 
Xinx 2 =x 

1 _ Ex lXae2 -mi(X,)ni 2 (X 2 ) 

x,nx 2 =0 



( 4 ) 



2.3. Dezert-Smarandache Theory (DSmT) 

In DSmT framework [34-36], the frame 0 = {0i, 0 2 , . . . , 0 n j is a finite set of n exhaustive elements which are not necessarily 
exclusive. The principle of the third excluded middle and Shafer’s model are refuted in DSmT (but can be introduced if needed 
depending on the model of the frame one wants to deal with). For a wide class of fusion problems, this is because the nature of 
hypotheses can be only a vague and imprecise or a crude approximation of the reality and no ultimate refinement is achievable. 
As a simple example, if we consider two suspects, Peter (P) and Mary (M), in some criminal investigations, it is possible that Peter 
has committed the crime alone, as well as Mary, or maybe Peter and Mary have committed the crime together. In that case, one 
has to consider the possibility forPn M ¥= 0, but there is no way to refine the original frame 0 = {P,M} into a finer one with exclu- 
sive finer elements, i.e. 0 = {P\(P n M),P n M,M\(P n M)}, because there is no physical meaning and no possible occurrence of 
the atomic granules P\(P n M) and M\(P n M). In other words, the exclusive elements of the refined frame satisfying Shafer’s 
model cannot always be well identified and precisely separated or may make no sense at all. This is why DSmT deals with non- 
exclusive, partially overlapped or vague elements and refutes Shafer’s model and third excluded middle assumptions. DSmT 
proposes to work in a fusion space defined by Dedekind’s lattice also called hyper-power set D®. 

The hyper-power set is defined as the set of all composite propositions built from elements of 0 with n and u operators 
such that [4]: 

(1) 0, fl,„. J n eD s ; 

(2) If A, Be D®, then AuBeD® andAnBeD 0 ; 

(3) No other elements belong to D®, except those obtained by using rules (1) or (2). 

Following Shafer’s idea, Dezert and Smarandache define a (generalized) basic belief assignment (or mass) as a mapping 
m( ): D® -> [0,1] such that: 

m(0) = 0 and y m(X) = 1 . 

XeD® 



Typically, when 0 = {A,B} and Shafer’s model doesnot hold, in DSmT one works with m(-) such that 

m(A) + m(B) + m(AuB) + m(AnB) = 1, (5) 

which appears actually as a direct and natural extension of (2) and (3). 

Actually DSmT also offers an advantage to work with Shafer’s model or with any hybrid model if some integrity con- 
straints between elements of the frame are known to be true and must be considered in the fusion. DSmT allows to solve 
static or dynamic 2 fusion problems in the same mathematical framework. For convenience, one denotes by G® the generalized 
fusion space or generalized power set including integrity constraints (that is, exclusivity as well as possible nonexistence 
restrictions between some elements of 0), so that G® = D® when no constraint enters in the model, or G® = 2® when one wants 
to work with Shafer’s model (see [34] for details and examples), or G® = 0 when working with probability model. If one wants 
to work with the space closed under union u, intersection n, and complementarity C operators, then G® = S®, i.e. the super- 
power set (see next section). A more general introduction of DSmT can be found in Chapter 1 of [36]. 

In DSmT, the fusion of two sources of evidence characterized by m i( ) and m 2 () is defined by m PC R 5 (0) = 0 and 
VXgG®\{0} 



tttpcj?5 (X) = m 12 (X) + ^ 



yeg® 

XnY=0 



' m,(X) 2 m 2 {Y) m 2 (X) 2 m t (Y) 

mi(X) + m 2 (Y) m 2 (X) + m,(Y) 



( 6 ) 



where all sets involved in the above formula are in the canonical form. m u (X) = m n (X ) = Z]x 1 x 2 ec® rn i(Xi)m 2 (X 2 ) corre- 

X t nX 2 =X 

sponds to the conjunctive consensus on X between n = 2 sources and where all denominators are different from zero. If a 
denominator is zero, that fraction is discarded. A general formula of PCR5 for the fusion of n > 2 sources has been proposed 
in [35]. 



i.e. when the frame or its model change with time. 
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2.4. Unification of fusion theory (UFT) 

Recently Smarandache has proposed in [33,35] an extension of DSmT by considering a super-power set S® as the Boolean 
algebra on 0, that is, S® = (6>,n,U,C(-)). In other words, S® is assumed to be closed under union u, intersection n, and com- 
plement C(-) of sets respectively. With respect to the partial ordering relation, the inclusion c , the minimum element is the 
empty set 0, and the maximal element is the total ignorance I = 0,-. Since it extends the power set space through the 

closed operation of n, U and complement (■) operators, that is, UFT not only considers the nonexclusive situation among 
the elements, but also considers the exclusive, exhaustive, nonexhaustive situation, and even the open and closed world. 
Typically, when 0 = {A,B}, in UFT one works with m(-) such that, 

m(A) + m(B) + m(A n B) + m(A u B) + m(C(A)) + m(C(B)) + m(C (A) u 0(B)) = 1. (7) 

3. Evidence Support Measure of Similarity (ESMS) 

Definition 3.1. Let’s consider a discrete and finite frame 0 and the fusion space C® including integrity constraints of the 
model associated with 0. The infinite set of basic belief assignments defined on G® is denoted by m G e. An Evidence Support 
Measure of Similarity (ESMS) of two (generalized) basic belief assignments rrq]-) and m 2 ( ) in m G e is the function 
Sim(-, ■) : m c e x m G e — > [0. 1] satisfying the following conditions: 

(1) Symmetry: Vmj(-),m 2 ( ■) € m C B,Sim(m^,m 2 ) = Sim(m 2 ,m i); 

(2) Consistency: Vm(-) e m G s,Sim(m,m) = 1; 

(3) Nonnegativity: Vrrii (•) , m 2 ( ■) 6 m G s,Sim(mi,m 2 ) > 0. 

This definition has been given in [16] and we focus on it from the point of view of belief function. We will say that m 2 ( ) is 
more similar to mi(-) than m 3 () if and only if Sim(m 1 ,m 2 ) > Sim(m^,m 3 ). The maximum degree of similarity is naturally ob- 
tained when both bba’s irq]-) and m 2 ( ) coincide, which is expressed by consistency condition (2). The equality 
Sim(m^,m 2 ) = 0 must be obtained when bba's of two different sources are not assigned to the same focal elements, that is, 
whenever irii(-) is focused on X 6 G®, which is denoted irif(-) and corresponds to m^X) = 1, and m 2 ( ) is focused on Y e G®, 
that is, m 2 (-) = mj(-) such that m 2 (Y) = 1, with Xn Y= 0. 

Obviously, for any bbami( ) e m c e (which is a |G®|-dimensional vector) and any small positive real number e, there exists 
at least one bba m 2 (-) 6 m G e for a given distance measure 3 d(-, ) such that d(rn 1 ,m 2 ) ^ e. 

Definition 3.2. (Agreement of evidence): If there exist two basic belief assignments m-i(-) and m 2 () in m c e such that for 
some distance measure d(-,-)> one has d(m 1 ,m 2 ) ^ e with e > 0, then nti(-) and m 2 {) are said e-consistent with respect to the 
distance d(-,-). 

As we know, the smaller e > 0 is, the closer the distance d(m 1 ,in 2 ) between m^) and m 2 () is, that is, the more similar or 
consistent nii(-) and m 2 ( ) are. 

From the previous definitions, ESMS is regarded as an interesting measure for evaluating the degree of similarity between 
two sources. We propose to use ESMS as a preprocessing/thresholding technique to reduce the complexity of the combina- 
tion of evidential sources by keeping in the fusion only the sources which are e-consistent. e is actually a threshold param- 
eter which has to be tuned by the system designer and which depends on the application and computational resources. 

4. Several possible ESMS 

In this section we propose several possibilities for choosing an ESMS function Sim(-, ) defined in the general fusion space 
G®, where all models, that is, probability, Dempster-Shafer (DS), Dezert-Smarandache (DSm) and Unification of Fusion may 
be chosen as the working model, however, the two evidential sources, when computing their similarity, must simultaneously 
work in the same model. The direct computation of their similarity between two sources respectively working in the differ- 
ent models is not significant. 

4.1. Euclidean ESMS function Sim E (m h m 2 ) 



Definition 4.1. Let 0 = {0^ 0 n ] (n > 1), rrq] ) and m 2 ( ) in m G s,Xj the ith (generic) element of G® and |G®| the cardinality of 

G®. The following Euclidean ESMS function can be defined 4 : 



3 Here we donot specify the distance measure and keep it only as a generic distance. Actually d(-, ) can be any distance measure. In practice, the Euclidean 
distance is often used. 

4 This is an extension of what is proposed in [13] for working in any fusion space, that is, in G 0 = 2 0 , G 0 = D 0 , or G & = S 0 . See [36], Chapter 1 for details. 
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Sim E (m u m 2 ) = 1 



x/2\ 



|C®| 

^(m 1 (X,-)-m 2 (X i )) 2 

i= 1 



( 8 ) 



Theorem 4.1. Sim E (m h m 2 ) defined in (8) is an ESMS function. 
(See proof in Appendix) 

4.2. Jousselme ESMS function Simj(m I ,m 2 ) 



Definition 4.2. Let m 1 ( ) and m 2 ( ) be two basic belief assignments in m c e provided by the sources of evidence St and S 2 . 
Assumed 5 a |G®| x \G @ \ positively definite matrix D = [Dy], where Dy = |X,- n X,-|/|X,- uXjj, with X,-, Xj e G®. Then, Jousselme ESMS 
function can be redefined from the Jousselme et al. measure [13]: 



Sirn^rn, , m 2 ) = 1 - ^(m, - m 2 ) T D(m, - m 2 ) 

or equivalently 

Simj(m i , m 2 ) = 1 - -J= + 2(m,,m 2 ), 

where (m 1 ,m 2 ) is the scalar product defined as: 



( 9 ) 



|G®[ |G®| 

(m,,m 2 ) = ^Dym,(X,)m 2 (Xj), 

i=i j=i 

X„ Xj g G®, i, j = 1 ,s, |G®|. 

Theorem 4.2. Simj(mt,m 2 ) defined in formula (9) is an ESMS function. 

(See proof in Appendix). 

Actually Sim E (m 1 ,m 2 ) is nothing but a special case of Sim / (m 1 ,m 2 ) when taking D as the |G®| x |G®|identity matrix. 
4.3. Ordered ESMS function Sim 0 (m,,m 2 ) 



The definition of the (partial) ordered ESMS function Sim 0 (m^,m 2 ) is similar to Simj(m^,m 2 ) but instead of using Jous- 
selme’s matrix D = [Dy], where Dy = |X,- n Xj|/|X,- u X/, with X,-, Xj e G®, we choose the DSm matrix 0= [Oy], where Oy = s(- 
X,- n Xj)/s(X, u Xj), and s(X) corresponds to the intrinsic informational content of the proposition X defined in details in [34] 
(Chapter 3). s(X) is used for partially ordering the elements of G®. More precisely, s(X) is the sum of the inverse of the length 
of the components of Smarandache’s code 6 of X. 

As a simple example, let’s take 0 = {0i,0 2 } with free DSm model (that is, when all elements are nonexclusive two-by-two), 
then the partially 7 ordered hyper-power set G® = D® is given by D® = {0,0 1 n 0 2 ,0i,0 2 ,0i u 0 2 }. As we know, s(0) = O, 
s(0i n 0 2 ) = 1/2, s(0i) = 1 + 1/2, s(0 2 ) = 1 + 1/2 and s(0j u 0 2 ) = 1 + 1 + 1/2 since Smarandache’s codes of 0, 8\, 0 2 , O-i n 0 2 and 
d-i u 8 2 are respectively given by {(■)} (empty code), {(1), (12)}, {(2), (12)}, {(12)} and {(1), (12), (2)}. The matrix 0 is defined by 8 



r s(0in« 2 ) 


s(«i n« 2 ) 


s(0,n0 2 ) 


s(0,n0 2 ) ] 












s(0,n« 2 ) 


s(f),) 


s(0 2 ) 


S(0,U0 2 ) 




■ 1 


1/3 


1/3 


1/5- 


s(0,no 2 ) 


s(0,) 


s(0,n0 2 ) 


s(«l) 




1/3 




1/5 


3/5 


s(«i) 


s(0,) 


S(0,U0 2 ) 


S(0,U0 2 ) 




1 


s(0,n« 2 ) 


s(0,n0 2 ) 


s ( 0 2) 


s(0 2 ) 




1/3 


1/5 


1 


3/5 


s (02) 


S(0! U0 2 ) 


S(0 2 ) 


S(0,U0 2 ) 






s(0,n« 2 ) 


s«h) 


s(0 2 ) 


s(0 1 uX 2 ) 




-1/5 


3/5 


3/5 


1 


|_S(0,Ufl 2 ) 


S(0,U0 2 ) 


S(0,U0 2 ) 


S(0,U0 2 ) J 













It is easy to verify that 0 is a positively definite matrix by checking the positivity of all its eigenvalues which are 
Aj = 0.800 > 0, X 2 rj 0.835 > 0, A 3 « 0.205 > 0 and X 4 as 2.160 > 0 respectively. We have verified that 0 is a positively definite 
matrix for Card(0) = n<5. However, it is an open challenge to prove that 0 is positively definite for any cardinality of 0 



5 Actually, Jousselme et al. in [13] did not prove that D = [D,j= X, o Xj /|X, u X/] is a positively definite matrix. D is until now assumed to be positively definite. 
This is only a conjecture and proving it is an open challenge. 

6 Smarandache code is a representation of disjoint parts of the Venn diagram of the frame 0 under consideration. This code depends on the model for 0. For 
example, let’s take 0 = {fli,0 2 }. If 0i n 0 2 = 0 (Shafer’s model) is assumed, then the code of 0i is (1), whereas if ()\ n (h / 0 (free DSm model) is assumed, then the 
code of 0i will be {(1), (12)}. The length of a component of a code is the number of characters between ( and ) in Smarandache’s notation. For example, the 
length of component (12) is 2. See [34], pp. 42-43 for details. 

7 This is a partial order since s(0 t ) = s(0 2 ). 

8 Actually, one works with C e \{0), and thus the column and row corresponding to the empty set do not enter in the definition of 0. 
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greater than 5, and that’s why we conjecture that 0 is a positively definite matrix for Card(0 ) > 5. Similarly to the formula 
(9), we define Szn7 0 (m 1 ,m 2 ) by: 

, 10 , 

It can be easily proved that Sim 0 (m^,m 2 ) is an ESMS function following a proof similar to the proof of Theorem 4.2. 

4.4. ESMS function Sim B (m h m 2 ) 

Another ESMS function based on Bhattacharya’s distance is defined as follows: 

Definition 4.3. Let mi( ), m 2 ( ) be two basic belief assignments in m c e, the ESMS function Sim B (m^,m 2 ) is defined by: 



Sim B (mi,m 2 ) = 1 



1 - J2 */m,(X,')m 2 (X,-), 



(ID 



where F is the core of sources Sr and S 2 , i.e. the set of elements of G® having a positive belief mass: F = {X e G®|mr(X) > 0 or 
m 2 (X) > 0}. 



Theorem 4.3. Sim B (m I ,m 2 ) defined in formula (11) is an ESMS function. 

(See proof in Appendix). 

5. Comparison of ESMS functions 

In this section we analyze the performances of the four ESMS functions aforementioned through a simple example, where 

0 = {0t,0 2 ,0 3 }. We assume the free DSm model for 0. In such case, G® = D® has 18 nonempty elements a,-, i = 1,2,. . ,,s 18. 

G® is closed under n and u operators according to Dedekind’s lattice. Of course, we also may choose other lattice structures 
in a similar way according to the right model of G®. 

Let’s assume that 0 2 is the true identity of the object under consideration. Its optimal belief assignment is denoted 
m 2 ( ) = {m 2 (0 2 ) = 1, m 2 (X) = 0 for X g C®\{0 2 }}. We perform a comparison of the four ESMS functions in order to show the 
evolution of the measure of similarity between mi( ) and m 2 () when m i(-) varies from an equally distributed bba to 
m 2 ( ). More precisely, we start our simulation by choosing m A (-) with all elements in D® equally distributed, i.e. 

m 1 (a,) = 1/18, for i=l,2,...,s, 18. Then, step-by-step we increase the mass of belief of 0 2 by a constant increment 

A = 0.01 until reaching m 1 (0 2 ) = l. Meanwhile, the masses m^X) of belief of all elements X^0 2 of G® take value 
[1 - m 1 (0 2 )]/17 to work with a normalized bba mr(-). The basic belief mass assigned to empty set is always zero, that is, 
m 1 (0) = m 2 (0) = 0. 

The degree of similarity of the four ESMS functions are plotted in Fig. 1. The speed of convergence 9 of a similarity measure 
is characterized by the slope angle a of the curve at origin, or by its tangent. Based on this speed of convergence criterion, the 
analysis of Fig. 1 yields the following remarks: 

(1) According to Fig. 1, tan(a B ) rs 0.86, tan(a £ ) rs 0.68, tan(a 0 ) ~ 0.6 and tan(o/) ss 0.57. Sim / (in 1 ,m 2 ) has the slowest con- 
vergence speed, then Sim 0 (m 1 ,m 2 ) takes second place. 

(2) Sim £ (mi,m 2 ) has a faster speed of convergence than Sim 0 (m 1 ,m 2 ) and Sim / (m 1 ,m 2 ) because it doesnot consider the 
intrinsic complexity of the elements in G®. 

(3) The speed of convergence of Sim B (m 1 ,m 2 ) is the fastest. When m i(-) and m 2 ( ) become very similar, Sim B (m 1 ,m 2 ) 
becomes very quickly close to 1. However, if a small dissimilarity between trii(-) and m 2 ( ) occurs, then Sim B (m^,m 2 ) 
becomes very small which actually makes it very sensitive to small dissimilarity perturbations. 

(4) In summary, one sees that no definitive conclusion about the best choice among these four ESMS functions can be 
drawn in general. However if one considers the speed of convergence as the evaluation criterion for the choice of a 
dissimilarity measure, one sees that Sim B (m 1 ,m 2 ) is the best choice, because it is very sensitive to small dissimilarity 
perturbations, whereas Simj(m 1 ,m 2 ) appears to be the worst choice with respect to such a criterion. 

6. Simulation results 

We present two simulation examples to show how ESMS filter performs within Generalized Fusion Machine (GFM), and 
to present its advantage. 



9 Here the convergence speed refers to how much the global agreement degree (similarity) is between m,( ) and m 2 () with the continuous decrease of 

mi(fl 2 ). 
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Fig. 1. Comparison of performance among ESMS functions. 



Table 1 

A list of given sources of evidence. 



s 


m(fli) 


m(e 2 ) 


m(e, n 6 2 ) 


m(0, u e 2 ) 


Sim E 


S, 


0.3 


0.4 


0.2 


0.1 


0.8735 


S 2 


0.3 


0.2 


0.4 


0.1 


0.800 


s 3 


0.4 


0.1 


0.2 


0.3 


0.8268 


s 4 


0.7 


0.1 


0.1 


0.1 


0.7592 


s 5 


0.1 


0.8 


0.1 


0.0 


0.5550 


s 6 


0.5 


0.2 


0.2 


0.1 


0.9106 


s 7 


0.4 


0.3 


0.1 


0.2 


0.9368 


S 8 


0.3 


0.1 


0.2 


0.4 


0.7592 


s 9 


0.4 


0.5 


0.1 


0.0 


0.8103 


Sio 


0.8 


0.1 


0.0 


0.1 


0.6806 


Sc, 


0.42 


0.28 


0.16 


0.14 


1.0000 


Sn 


0.5 


0.0 


0.2 


0.3 


0.7569 


$12 


0.2 


0.6 


0.1 


0.1 


0.7205 


$13 


0.4 


0.3 


0.2 


0.1 


0.9360 


S14 


0.9 


0.1 


0.0 


0.0 


0.6230 


Si 5 


0.5 


0.2 


0.1 


0.2 


0.9100 


Si 6 


0.5 


0.3 


0.0 


0.2 


0.8900 


Si 7 


0.5 


0.0 


0.1 


0.4 


0.7205 


Si 8 


0.7 


0.2 


0.1 


0.0 


0.7807 


Si 9 


0.1 


0.7 


0.1 


0.1 


0.6217 


S 20 


0.3 


0.6 


0.1 


0.0 


0.7390 


Sc2 


0.44 


0.29 


0.13 


0.14 


1.0000 



Example 1 . Let’s take a 2D frame of discernment 0 = {Oi,0 2 } and consider 20 equireliable sources of evidence according to 
Table 1. We consider the free DSm model and the fusion space is the hyper-power set D® = { 61 , 62,61 D d 2 , u 0 2 }. S cl denotes 
the barycentre of the front ten belief masses, while S c2 denotes the barycentre 10 of all belief masses. The measure of similarity 
based on Euclidean ESMS function defined in Eq. (8) has been used here, but any other measures of similarity could be used 
instead. In this example, if we take 0.75 as the threshold value, we see from Table 1 and for the former 10 sources of evidence, 
that the measures of similarity of S 5 and S 10 with respect to S cl are lower than 0.75. Therefore, the sources S 5 and S 10 will be 
discarded/filtered in the fusion. If the threshold value is set to 0.8, then the sources S 5 , S ]0 , S 4 and S 8 will be discarded. That is, 
the higher the given threshold is, the fewer the number of information sources through the filter is. 

The main steps of the algorithm for carrying out GFM are the following ones: 

(1) Initialization of the parameters: the number of sources of evidence is set to zero, (that is, one has initially no source, 
s = 0), so that the number of sources in the filter window is n = 0. 



10 Let’s denote k = |C®| the cardinality of G® and consider S independent sources of evidence. If all sources are equireliable, the barycentre of belief masses of 
the S sources is given by: Vj = 1 k, m(Xj) = see [21 ] for details. 
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Fig. 2. Fusion result of the front 10 sources using different thresholds (0 0.9). 

(2) Include 11 a source of evidence S s and then test if the number of sources s is fewer than two. If s > 2, then go to the next 
step, otherwise include/take into account another source of evidence S 2 . 

(3) Based on the barycentre of gbba’s of the front n ^ 10 evidence sources, the degrees of similarity are computed accord- 
ing to the formula (8), and compared with a prior tuned threshold. If it is larger than the threshold, then let n = n + 1. 
Otherwise, introduce a new source of evidence S s+1 . 

(4) If n = 1, the current source, say S, is not involved in the fusion. If n = 2, then the classical DSm rule [34] is applied in the 
fusion step between S and S 2 , that is, the conjunctive consensus. Then PCR5 rule [35] is applied to redistribute the 
remaining partial conflicts only to the sets involved in the corresponding partial conflicts. We get a new combined 
source affected with the same index S. If 2 < n ^ 10, after the current evidence source S s is combined with the final 
source of evidence produced last time, a new source of evidence is obtained and assigned to S again. Whenever 
n ^ 10, go back to step (2), otherwise, the current source of evidence S s under test has been accepted by the ESMS fil- 
ter, Si is assigned to S,-_ i, i e [2,s, 10], and S s is assigned to S] 0 . Then, Si 0 is combined with the last source S, the com- 
bined result is reassigned 12 to S, and then, go back to step (2). 

(5) Test whether to stop or not 13 : if not, then introduce a new source of evidence S s+ i, otherwise stop and exit. 

We show two simulation results in Figs. 2 and 3 following the working principle of GFM, when we use the sources of 
evidence listed in Table 1. 

The comparison of Figs. 2 and 3 yields the following remarks: 

(1 ) On Fig. 2, we donot see the real advantage of ESMS filter since the convergence to di without ESMS filter (green curve) 
is better than with ESMS filter in terms of improving fusion precision. This is because some useful sources of informa- 
tion are filtered and thrown away with the increase of the threshold. Therefore the number of sources entering in the 
final fusion will decrease. This yields a slower speed of convergence to ffi. For example, let’s consider the source S 3 in 
Table 1. If S 3 is combined with itself only once, according to formula (6), the fusion result is S = [m^ffi) = 0.5707, 
m N (0 2 ) = 0.0993, m N (0i n 0 2 ) = 0.1680, m N (0i u 0 2 ) = 0.1620]. 14 If combined twice, the result is S= [0.6765, 
0.0852,0.1429,0.0954]. If combined thrice, then the result is S = [0.7453,0.00693,0.1216,0.0638]. The more the combina- 
tional times is, the nearer to 1 m N (0x) is, and the nearer to 0 m N {0 2 ) is. Therefore, ESMS filter might also result in losing 
some useful information, while it filters some bad information. 

(2) On Fig. 3, one sees the role played by ESMS filter. When there are highly conflicting sources, the result of the fusion 
will not converge if ESMS filter is not used. With the fine-tuning of ESMS threshold, the convergence becomes better 
and better because ESMS filter processes the fused information, and withdraws the sources which might cause the 
results to be incorrect or imprecise, so that it improves the fusion precision and correctness. 



11 We assume the free DSm model and consider that the general basic belief assignments are given. 

12 In this work, we also use an ESMS filter window in a sliding mode. 

13 in our experiment, judge whether the mobile robot stops receiving sonar’s data. 

14 m N ( ) refers to the new generalized basic belief assignment. 
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Fig. 3. Fusion result of the total 20 sources using different thresholds (0 ~ 0.9). 



(3) We increase the applicability of the classical rules of combination. Since we reduce the number of conflicting sources 
of evidence combined thanks to ESMS preprocessing, the degree of conflict between sources combined is kept low. 
Therefore, classical rules, like Dempster’s rule, which do not perform well in high conflicting situations can also be 
applied in robot perception. Without ESMS preprocessing step, the classical fusion rules cannot work very well 
[26,32]. Therefore, we extend their domain of applicability when using ESMS filtering step. 



Example 2. To show the advantage of the proposed method, we compare the performance of the method 1 (that is, using 
sequential fusion with PCR5 only) with the method 2 (using sequential fusion PCR5 with ESMS preprocessing) in terms of 
complexity cut. We consider the following two different cases: 

• Case 1 : (simple) 

We consider the frame 0 = {0i,0 2 } with DSm model. There are four focal elements in D®\{0}. We will consider 1001 
sources of evidence having the following same belief assignments: m s (0i) = O.4, m s (0 2 ) = 0.2, m s (0i n 0 2 ) = 0.3, 
m s (0-, u 0 2 ) = 0.1, 1 ^ s ^ 1001. 

• Case 2: (a bit more complex) 

We consider a bit more complex case with 0 = {0!,0 2 ,0 3 } with DSm model. There are 18 elements in D®\{0). If all ele- 
ments of D®\{0} are focal elements, then the computation of the combination becomes very costly and very time-con- 
suming. We only consider 1001 sources of evidence having the following same belief assignments over partial 
elements in D®\{0}: m s (0i) = 0.15, m s (0 2 ) = 0.25, m s (0 3 ) = O.l, 015(0! n 0 2 ) = 0.05, rn s (0i n 0 3 ) = 0. 08, m s (0 2 n 0 3 ) = 0.18, 
m s (0 1 n 0 2 n 0 3 ) = 0.12, m s (0 1 u 0 2 u 0 3 ) = 0.07, 1 ^ s ^ 1001. 

In either Case 1 or 2, all sources are the same. Of course, we had better produce a set of sources (1001) at random, and 
then partially choose sources by setting different ESMS thresholds. But a set of different sources (1001) produced at random 
are not easy to be listed here. Since all operations involved in each method do count in the computation time, that is, 
computation of the barycenter of belief masses, computation of ESMS, comparison and combination operations, obviously, 
there is nearly no difference between the 1001 same and different sources in term of computational burden. For 
convenience, here we replace the 1001 different sources with the 1001 same ones. 

Table 2 shows, for Case 1, the comparison of the computing time 15 for the sequential fusion based on PCR5 only 16 (method 
1) with the sequential fusion using PCR5 and ESMS preprocessor (method 2) when considering 21, 51 101, 201, 401, 501, 801 
and 1001 sources of evidence selected. Similarly, Table 3 shows the comparison of computing time for Case 2. 

The comparison of Tables 2 and 3 yields the following remarks: 



15 using Processor: Intel (R) Core (TM)2 Duo CPU E7500 @ 2.93 GHz. EMS memory: 2.0 GB. 

16 That is without ESMS preprocessing. 
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Table 2 

Computing time with methods 1 and 2 for Case 1 . 



Number of sources 


1001 


21 


51 


101 


201 


401 


501 


801 


1001 


Time in ps 


241 


167 


172 


181 


199 


233 


248 


300 


335 


Method used 


1 


2 


2 


2 


2 


2 


2 


2 


2 



Table 3 

Computing time with methods 1 and 2 for Case 2. 



Number of sources 


1001 


21 


51 


101 


201 


401 


501 


801 


1001 


Time in ps 


916 


231 


255 


295 


376 


536 


619 


864 


1025 


Method used 


1 


2 


2 


2 


2 


2 


2 


2 


2 



(1) Even if we have introduced an ESMS preprocessing step (method 2), it turns out that finally a drastic cut of computing 
burden is obtained when these sources of evidence are highly conflicting and inconsistent. This is because we can sig- 
nificantly reduce the number of sources combined with ESMS criterion. Seen from Table 2 or Table 3, when 21 sources 
are chosen from 1001, that is, the remaining 800 sources assumed/regarded as highly conflicting and inconsistent 
sources (actually they are consistent absolutely, here it is only a hypothesis), the computing time is 167 ps in Table 2 
and 231 ps in Table 3, whereas the computing time using method 1 is 241 ps in Table 2 and 916 ps in Table 3. There- 
fore, we save 1 - (167/241 ) « 30% of computing time in Table 2 when combining only 21 sources among 1001, 
whereas we save 1 - (231/916) s» 75% of computing time in Table 3. 

(2) If the sources of evidence are less conflicting and more consistent, then ESMS preprocessing step seems useless. 
Indeed, one sees in Table 2 that the computing time using method 2 is 335 ps after combining 1001 sources, which 
is greater than the computing time of the method 1. The same remark holds for Case 2 as shown in Table 3. In such 
a very particular case (when no source is filtered through ESMS), the computing time with method 2 becomes greater 
than that of method 1 , because method 2 needs more computation to carry out ESMS preprocessing without gaining 
benefit (that is, saving computations) in PCR5 fusion step. Of course, with the increase of the number of focal elements, 
the computations for PCR5 fusion become more and more complex and the time taken by ESMS preprocessor becomes 
negligible with respect to PCR5 fusion step. 



7. An application in mobile robot perception 



The information acquired for building grid map using sonar sensors on a mobile robot is usually uncertain, imprecise and 
even highly conflicting. An application in autonomous robot perception and navigation provides a good platform to verify 
experimentally the advantage of ESMS filter in GFM. Although there exist many methods of building map based either on 
Probability Theory [39], FST (Fuzzy System Theory) [27], DST [37], GST (Grey System Theory) [42-44], or DSmT [23], we just 
compare the performances of the map building using a classical fusion machine without ESMS filter (i.e. CFMW) with the 
classical fusion machine with ESMS filter (called GFM) in DSmT framework only. A detailed comparison between our current 
ESMS-based approach with other methods is given in a companion paper in [24] where we show that ESMS-based approach 
outperforms other approaches using almost the same experimental conditions and inputs. In order to further lessen the mea- 
surement noises, we improve our past belief assignment model of sonar sensors in DSmT framework 17 as follows: 



m(9 1 ) = 



m{0 2 ) = 



(1 ~ P/{R~2t)) x (1 - 2/2) if {*' 

0 otherwise, 



Rmin ^ P ^ R t . 
^ (p ^ CO/ 2, 



exp(-3(p - R) 2 ) x 2 if 

0 otherwise, 



Rmin ^ ^ H T 6, 

0 ^ ^ co/2, 



m(0, n B 2 ) = 



l-(2 ( P -(R-e)/R) 2 ) if 

0 otherwise. 



Rmin ^ R + €, 

0 < cp < 00/ 2, 



( 12 ) 



(13) 



(14) 



17 We assume that there are only two focal elements 9 i and 0 2 in the frame of discernment. Elements of hyper-power set are 0\, 6 2 , 9 \ n 0 2 and 9 \ u 0 2 . 9 \ 
represents the emptiness of a given grid cell, 0 2 represents the occupancy for a given grid cell, 9 ^ n 9 2 means that there is some conflict between two sonar 
measurements for the same grid cell and 9 \ u 9 2 represents the ignorance for a grid cell because of the possible lack of measurement. 
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Fig. 4. The running robot in a real experimental environment. 



m{9 1 
where, 2 is 



u02) = ( tan(2 (p-R)) x (1 -2) 

\ 0 otherwise, 

given by (see [8] for justification) 



if 



ft < p < ft + e, 



2 = 



l-(2<p/co) 2 if 0<|<p|< co/2, 
0 otherwise. 



(15) 



(16) 



The parameters ft, p, e, ft m , n , co, and (p in formulas (12)-(16) were defined and used in [22-24], ft is the sonar reading (mea- 
surement). p is the distance between the grid cell and sonar’s emitting point, e is the measurement error. ft m!n is the minimal 
range of sonar sensors. <±> is the scattering angle of sonar. <p is the angle between the line (from the grid cell to sonar emitting 
point) and the sonar’s emitting direction. The following functions Cj and C 2 play an important role in lessening noises in the 
process of map building. (j function, proposed by Wang in [42-44], is a constriction function 18 for sonar measurements de- 
fined by: 



C, 



'0 

Pi,-p 

' Ph~Ph 

1 



if P>Pi 2 , 
if Pi,<P<Pi 2 , 
if P<Pi 2 , 



(17) 



where, p , and p, 2 represents the upper and lower limits of valid measurements, C 2 is the constraint function 19 for sonar’s 
uncertainty defined as follows: 



C 2 



(P-R- i-0.5e \2 

V 0.5e ) 

(p-R- 0.5e \2 

V 0.5e / 

o 



if P ~ 

if P ~ 
if \P — 



ft> -0.5e, 
ft < 0.5e, 
ft| > 0.5e, 



(18) 



where, the product of Cj and C 2 is multiplied by the belief assignment function (m(0 2 ),m(0 1 n 0 2 ),m(0 1 u 0 2 )) respectively. 

The experiment is performed by running a Pioneer II mobile robot with 16 sonar detectors in the indoor laboratory envi- 
ronment as shown in Fig. 4. The environment’s size is 4550 mm x 3750 mm. The environment is divided into 91 x 75 rect- 
angular cells having the same size according to the grid map method. The robot starts to move from the location (1 m,0.6 m), 
which faces towards 0 degrees. We take the left bottom corner as the global coordinate origin of the map. Objects/obstacles 
in the rectangular grid map are shown in Fig. 5. The processing steps of our intelligent perception and fusion system have 
been carried out with our software Toolbox developed under VC++ 6.0 and with OpenGL server as a client end. When the 
robot moves in the environment, the server end collects much information (such as the location of robot, sensors measure- 
ments, velocity) from the mobile robot and its sensors onboard. Through the protocol of TCP/IP, the client end can get any 
information from the server end and fuse them before displaying final result (the estimated grid map). 



18 The main idea is that the sonar readings must be discounted according to sonar characteristics. 

19 The main idea consists in assigning high belief assignments to sonar readings close to the sonar sensor. 
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mobile robot 
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1 000mm 



Environment 



1 100mm 



1800mm 



3750mm 

Fig. 5. Global coordinate system for the experiment. 



Y 




3750mm 



Fig. 6. Map building based on GFM before improving the sonar model. 



Since our environment is small, the robot moves on a short distance during a relatively short period of time. Therefore, 
one only considers the self-localization method based on <5-NFAM 20 method [20,22] with the search from 0 - S 0 to 0 + S 0 . In 



20 One has taken S 0 = 5° in our experiment. 
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Fig. 7. Map building based on GFM after improving the sonar model. 



Y 




3750mm 



Fig. 8. Map building based on CFMW before improving the sonar model. 
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order to reduce the computation burden, the restricted spreading arithmetic has been used. The main steps of map building 
with GFM are the following ones: 

(1) Initialize the parameters of the robot (location, velocity, etc.). 

(2) Acquire 16 sonar measurements, and robot’s location from the odometer, when the robot is moving in the environ- 
ment. The robot’s pose is calibrated with our <5-NFAM method [20,22] and the time clock with a 100 ms period is 
set up. 

(3) Compute gbba of the fan-form area detected by each sonar sensor according to the formulas in [23]. 

(4) Apply DSmT-based GFM, that is, adopt Euclidean information filter to choose basic consistent sources of evidence 
according to the formula (8). Then combine the consistent sources with DSm conjunctive rule [3,4,34] and compute 
gbbas after combination. Then, redistribute partial conflicting masses to the gbba’s of sets involved in the partial con- 
flict only with PCR5 rule [35]. 

(5) Compute the belief of occupancy Bel(0 2 ) of some grid cells according to [34]. Save them into the map matrix and then 
go to step (6). 

(6) Update the map of the environment (here we set the second timer, of whose interval is 100 ms). Generally speaking, 
the more the times of scanning map are, the more accurate the final map rebuilt is. At the same time, also test whether 
the robot stops receiving the sensed data: if yes, then stop fusion and exit, otherwise, go back to step (2). 

In this experiment, we obtain the maps built by GFM before and after improving the sonar model as shown in Fig. 6 and 7 
respectively. In order to show the advantage of ESMS filter in GFM, we also compare our approach with the classical fusion 
machine without ESMS filter (called CFMW). The maps built by CFMW before and after improving the sonar model are 
shown in Figs. 8 and 9 respectively. Whenever the map is built before or after improving the sonar model, one sees that 
GFM always outperforms CFMW because one obtains clearer boundary outlines and fewer noises in the map built, in addi- 
tion, ESMS information filter coupled with PCR5 fusion rule, allows to reduce drastically the computational burden because 
ESMS filter can filter the outlier-sources. With GFM approach, only the most consistent sources of evidence are combined 
and this allows to reduce the uncertainty in the fusion result and to improve the robot perception of the surrounded 
environment. 



Y 




3750mm 



Fig. 9. Map building based on CFMW after improving the sonar model. 
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8. Conclusions 

In this paper, a general evidence supporting measure of similarity (ESMS) between two basic belief assignments has been 
proposed. ESMS can be used in different fusion spaces (lattice structures) and with different distance measures. This ap- 
proach allows to select the most coherent subset of sources of evidence available and to reject outlier-sources which is seri- 
ously inconsistent with other sources. Therefore, a drastic cut of computational burden is possible with keeping good 
performances which is very attractive for real-time applications having limited computing resources. The hybrid of ESMS 
with the sophisticated and efficient PCR5 fusion rule of DSmT, called GFM (Generalized Fusion Machine), is specially useful 
and interesting in robotic applications involving real-time perception and navigation systems. The real application of GFM 
for mobile robot perception from sonar sensors presented in this work shows clearly a substantial improvement of the fusion 
result in map building/estimation of the surrounded environment. This work also shows an important role played by the 
most advanced fusion techniques for applications in robotics. 
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Appendix A 



Theorem 4.1. Sim E (m h m 2 ) defined in (8) is an ESMS function. 

Proof 

(1) Let’s prove that Sim E (m^,m 2 ) e [0, 1 ].If Sim £ (m 1 ,m 2 ) > 1, from (8) one would get — m 2 (X,-)) 2 < 0 which 

is impossible, so that Sim £ (m 1 ,m 2 ) < 1. Let’s prove Sim £ (m 1 ,m 2 ) > 0 or equivalently from (8), Jj|f^(mi(X,) — m 2 (X ,)) 2 ^ 
2.This inequality is equivalent to 1 mi (X,) 2 + S£i l m 2 (X,) 2 C 2 + 2)T)J f '” l m l (X,j m 2 (X,).We denote it (i) for short.(i) 
always holds because one has (£Ei lm i(Xi) 2 + £i=i lm 2 (Xi) 2 ) < ( [Xjl=i lm i(X,)] + [eEi 1 m 2 (X,-) j ^ and thus 

(Ew'miCXj) 2 + Ei=i l m 2 (X i ) 2 ) ^2 because ^] r;f | <l m s (X i )| =1 for s=l,2 (m s () being normalized bba).Therefore 
inequality (i) holds and thusSim £ (m 1 ,m 2 ) ^ 0. 

(2) It is easy to check that Sim £ (m 1 ,m 2 ) satisfies the first condition of Definition 3.1. 

(3) If mfi) = m 2 (), then Sim £ (m 1 ,m 2 ) = 1 because 

|C®| 

£(m,(X0-m 2 (X,)) 2 = O. 

i=i 

Thus the second condition of Definition 3.1 is also satisfied. 

(4) Nonnegativity has been proven above in the first part. Herein we use a particular case to show thatSim(ml,m2) = 0, i.e. 
there exist mf and mj for some X, Y e G <s, \{0} such that X^Y, then according to (8), one gets 
£)|f^(mi(Xj) - m 2 (X,)) 2 = [mf(X)] 2 + [mJ(Y)] 2 = 2 and thus one has Sim £ (mf, m%) = 1 - (\/2/\/2) = 0, so that Sim £ (-, ) 
verifies the third condition of Definition 3.1. □ 



Theorem 4.2. Simj(m I ,m 2 ) defined in formula (9) is an ESMS function. 

Proof 

(1) Since the matrix D is conjectured to be a positively definite matrix, Simj(m^,m 2 ) satisfies the condition of symmetry. 

(2) If mi is equal to m 2 , according to (9), one getsSim;(mi,m 2 ) = 1. On the other hand, if Sim/(m 1 ,m 2 ) = 1, then the condition 
mi = m 2 holds. That is, the condition of consistency is satisfied. 

(3) According to (9), it can be drawn that Sim/(m 1 ,m 2 ) ^ Sim £ (m 1 ,m 2 ), and since the minimum value of Simj(m^,m 2 ) is zero, 
then Sim J (m 1 ,m 2 ) is non-negative. 
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(4) According to the definition of Simj(m 1 ,m 2 ), we can easily verify that Simj(m^,m 2 ) is a true distance measure between 
m, and m 2 . □ 

Theorem 4.3. Sim B (m I ,m 2 ) defined in formula ('ll) is an ESMS function. 

Proof 

(1) Since J2xerV m '(Xi)m 2 (Xi) = Sx-eF \/ m 2 (-Xf)mi(X;), then Sim B (m 1 ,m 2 ) satisfies the condition of symmetry. 

(2) If mi(-) = m 2 ( ), according to (11), 

22 v^i (Xi)m 2 (Xi) = 22 m i( x 0 = 1 

X,eF X,£F 

and therefore Sim B ( mi, mi) = 1. On the other hand, if Sim B (m 1 ,m 2 ) = 1, then the condition m i(-) = m 2 ( ) holds. That is, the 
condition of consistency is satisfied. 

(3) From the definition of bba, F m i(X) = 1- Therefore, 

22 \Zmi(Xi)m 2 (Xi) e [0, 1], 

X,£F 

According to (11), it can be drawn that Sim B (m 1 ,m 2 ) e [0,1]; that is, the minimum value of Sim B (m 1 ,m 2 ) is zero. There- 
fore, Sim B (m 1 ,m 2 ) is a nonnegative measure. 

(4) According to the definition of Sim B (m 1 ,m 2 ), we can easily verify that Sim B (m 1 ,m 2 ) is a true distance measure between 

mi and m 2 . □ 
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